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Due to practical reasons, experimental and theoretical continuous-variable (CV) quantum information (QI) 
has been heavily based on Gaussian states. Nevertheless, many CV-QI protocols require the use of non-Gaussian 
states and operations. Here, we show that the Einstein-Podolsky-Rosen steering inequality can be used to obtain 
a practical witness for the generation of pure bipartite non-Gaussian states. While the scenario require pure 
states, we show its broad relevance by reporting the experimental observation of the non-Gaussianity of the CV 
two-photon state generated in the process of spontaneous parametric down-conversion (SPDC). The observed 
non-Gaussianity is due only to the intrinsic phase-matching conditions of SPDC. 

PACS numbers: 42.50.Dv 


I. INTRODUCTION 


II. THE EPR-STEERING INEQUALITY WITH BIPARTITE 
GAUSSIAN STATES 


Continuous-variable (CV) quantum information (QI) is a 
research field that has been increasingly growing in the past 
few years [1], The need to cover larger Hilbert spaces is 
motivated by QI protocols, such as quantum key distribution 
(QKD), that has advantages when implemented in higher di¬ 
mensions [2, 3]. The extension of many QI protocols, first 
proposed considering discrete quantum systems, to the realm 
of CV systems has been heavily based on Gaussian states [4— 
6]. This is due to the fact that their covariance matrices are 
fully determined by the first and second order moments [7, 8], 
and also because of the practicality in the creation, manipula¬ 
tion and detection of Gaussian states [9, 10]. 

Nevertheless, several recent works have shown the rele¬ 
vance of non-Gaussian states and operations [11]. For in¬ 
stance, they are required for quantum computation with clus¬ 
ter states [12], entanglement distillation [13, 14], quantum er¬ 
ror correction [15] and loophole-free Bell tests [16, 17]. Be¬ 
sides, they provide advantages to the quantum teleportation, 
quantum cloning and state estimation tasks [19-22]. 

In this work, we show that the Einstein-Podolsky-Rosen 
(EPR) steering inequality [23-26] can be used to obtain a 
witness for the non-Gaussianity of pure bipartite CV quan¬ 
tum states. While the scenario require pure states, we show 
its broad relevance by reporting the observation of the non- 
Gaussianity of the CV two-photon state generated in the pro¬ 
cess of spontaneous parametric down-conversion (SPDC)[27- 
30]. SPDC is up to date the most used source for experimen¬ 
tal investigations in the field of quantum information, and our 
work highlights the simplicity of using this source for new ap¬ 
plications in CV-QI. The generated down-converted photons 
are correlated in their transverse momenta and can be used 
to test the EPR-paradox [31]. The observed non-Gaussianity 
is due only to the intrinsic phase-matching conditions of the 
SPDC process [32, 33], thus, highlighting the simplicity of 
using SPDC sources for new applications of CV-QI. 
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Consider a bipartite system described by a pure CV-state 
P 12 = \v)(¥l ant l a generic pair of complementary non¬ 
commuting observables if and v,-, with i = 1,2 used to de¬ 
note the operation at each subsystem. The spectral decom¬ 
position of if and v, is a infinite set of continuous variables. 
In Ref. [25] it was introduced the EPR-steering criterion as 
Aj n f{u 2 )tf n f(v 2 ) > C, where C is is a value that depends on 
the chosen observables. The violation of this inequality im¬ 
plies the implementation of a EPR-paradox [31, 34]. The in¬ 
ferred variances are given by 

A ?»/(“ 2) = J duiP(ui)A 2 (u 2 \ui), (1) 


where A 2 (u 2 \l<i) is the variance of the conditional probability 
distribution P(u 2 \u\), and P(u\) is the marginal probability 
distribution of one party’s outcomes. 

Now, let us consider a general continuous variable 
two mode pure state described by a Gaussian amplitude 
/// (q |. q2 ), where the vectors qi and q 2 are the outputs of the 
observable q on each party 


^ b(qi^ 2 ) exp 
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exp 
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( 2 ) 


with <7+ and <7 being the widths of the corresponding Gaus¬ 
sian functions. Several CV physical systems can be modelled 
in this way, for example: an atomic ensemble interacting with 
an electromagnetic held, two entangled photons sent through 
Gaussian channels, photoionization of atoms and photodisso¬ 
ciation of molecules [35], spontaneous emission of a photon 
by an atom [36-38], and multiphoton pair production [39,40]. 
Since Gaussian amplitudes have the same functional form in 
both transverse directions we can, without loss of generality, 
work in one dimension and consider only its scalar form. 

From Eq. (1) one can obtain the limit of C for Eq. (2) 
[40], Note that A 2 nfG (q 2 ) = A 2 (q 2 \q { ) = a 2 C7 2 /(°+ + <*-)■ 
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One can also find the inferred variance of the complemen¬ 
tary observable and it is given by A 2 , c (x 2 ) = A 2 (x 2 |xi) = 
1 / (<T+. + C 2 ). Therefore, the EPR-steering inequality reads 

A 2 (?2ki)A 2 (x 2 |xi) = -^±—_ > I. ( 3 ) 

(a 2 +a 2 ) 2 4 

Let us define one parameter P = <t_/<t for simplicity. The 
Schmidt Number Kq for the general Gaussian state of Eq. (2) 
is given by [41] 

k ° =l S +p ) 2 ■ (4) 

Thus, the EPR-steering inequality can be written as 

A 2 {q 2 \qi)A 2 {x 2 \xi) = J (5) 

Note that due to the symmetry of Eq. (2) one has that Kq = 
KQ a x K G i 3 , where Kqj represents the Schmidt number in the 
transverse direction j = a, fi [42]. Thus, one may further sim¬ 
plify the EPR-steering inequality for Gaussian states to 

W = A 2 (q 2 \qi)A 2 (x 2 \xi) = J (6) 

Got 

Thus, one may test it by performing the measurements in only 
one transverse direction. Moreover, it does not depend of 
the values chosen for q\ and jq. When P = 1, Kq = Kq u = 
Kcp = 1 and the Gaussian state is a product state. In Fig. 1(a) 
[Fig. 1(b)] we show Kq u ( W ) with a solid blue line, while 
varying P. 


III. A WITNESS FOR THE NON-GAUSSIANITY OF CV 
QUANTUM STATES 

From the results obtained above it is possible to envisage 
a simple and practical way to determine if a certain bipar¬ 
tite pure state is Gaussian or not. Note from Eq. ( 6 ) and 
Fig. 1(b) that pure bipartite Gaussian states will always vio¬ 
late the EPR-steering inequality. The reason is that they are in 
general entangled states [See Fig. 1(a)]. The only exception is 
the point marked with the horizontal dashed line in Fig. 1(b), 
which represents the point where the Gaussian state is a prod¬ 
uct state. In this case, we have that W = = 0.25, which 

corresponds to the upper quantum bound for pure bipartite 
Gaussian states. Thus, the observation of a value greater than 
0.25 for W with pure entangled states can only be achieved 
while considering non-Gaussian states. 

While the scenario require pure states, now we show its 
broad relevance by reporting the observation of the non- 
Gaussianity of the CV spatial two-photon state generated in 
the process of SPDC [27-30]. When perfect colinear phase¬ 
matching is considered and neglecting effects of anysotropy, 
we can write the spatial two-photon state as [28, 29] 

IV / ')i2 oc Jjdqidq 2 E P (qi + q 2 )G(qi —q 2 )|lqi)|lq2), (7) 


(a) 




(b) 



P 



FIG. 1. (Color Online) (a) Shows the Schmidt number for bipartite 
Gaussian states and for the CV spatial state of SPDC [Eq. (7)], while 
varying P and considering one transverse direction of q^. (b) The 
values of W plotted in terms of P. 


where |lq.‘> represents one photon in mode i (i = 1 , 2 ) usu¬ 
ally called signal or idler, and with the transverse momen¬ 
tum q. E p ( q) is the angular spectrum of the pump beam. 
Usual experimental configurations adopt a Gaussian pump 
beam and in this case E p ( q) « exp [—c 2 |q| 2 /4]. c represents 
the beam radius at the crystal plane. G(q) = sine (Z?|q| 2 ) 
defines the phase-matching conditions of the SPDC process, 
with sinc(^) = Sln |^ . b is defined by b = ^ , where L is the 
crystal length, and k the wavenumber of the down-converted 
photons. In terms of these definitions, P reads i \J~^- This 
state can be rewritten in the complementary transverse posi¬ 
tion representation as [33, 43] 

Ir)i2« JJdx\dxi E p ^ J X2 ^j G ^ 2 * 2 ^ |lxi)|l* 2 ), 

( 8 ) 

where the functions E p (x) and G(x) are the Fourier transform 
of E p { q) and G(q). Thus, E p (x) exp [—|xr| 2 /c 2 ] and G(x) 

1 - i Si (45 M 2 ) = sint (4b |*| 2 ), where Si(x) = /q dt sinc(f). 
Here, the transverse position xj and transverse momentum qj 
are the complementary observables for the EPR-steering in¬ 
equality test [31]. Clearly, the CV spatial two-photon state of 
SPDC is a non-Gaussian state [32]. 

An experimental observation of the effects that arise form 
the phase-matching conditions, namely, the non-Gaussianity 
of |i//)i 2 , was reported in Ref. [33]. They demonstrated how 
the spatial correlations in the near field plane of a non-linear 
crystal changes when the phase matching conditions varies. 
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FIG. 2. (Color Online) Experimental setup, (a) State preparation stage, (b) Setup configuration for measuring A 2 (^ 2 1^1) ■ ( c ) Setup configura¬ 
tion for measuring A 2 (x 2 \xi). See the main text for details. 


Now, we demonstrate experimentally how W can be used to 
detect the non-Gaussianity of this state. The demonstration 
is based on the fact that the spatial state of SPDC process is 
pure and entangled [41], even considering a post-selected one 
transverse direction [30, 32], This means that for any value 
of P, the Schmidt number K S a is always greater than 1. In 
Fig. 1(a) [Fig. 1(b)] we show K$ a ( W ) with a dashed red line, 
while varying P for the state |v / )i 2 - One can see that for some 
values of P (0.56 < P < 2.58) the values of W are greater than 
3 , thus, witnessing the non-Gaussianity of this state. 


down-converted photons at 710 nm with 5 run of bandwidth. 
This fact implies no entanglement between the frequencies 
and the transverse spatial coordinates, and then the reduced 
spatial state must be pure. Besides, the usage of spatial filters 
is explained afterwards. 

Furthermore, a polarizer beam splitter (PBS) separates the 
signal and idler photon modes. For measuring A 2 (^ 2 1^ 1 ) and 
A 2 (a' 2 |a'i ), we performed conditional coincidence measure¬ 
ments between the idler and signal photons in two different 
transverse planes (See Fig. 2): the far- and near-field planes 
of the non-linear crystal, respectively [31]. 


IV. EXPERIMENT 


A. Measurement of A 2 (c/ 2 \q \) 


The experimental setup is illustrated in Fig. 2. We used a 
solid-state laser source, at 355 nm, to pump a /3-barium-borate 
type-II (BBO-II) non-linear crystal (L = 1.8 cm) for the gen¬ 
eration of the down-converted photons. Initially, the Gaus¬ 
sian pump beam had a waist of c = 200 / 1 /// at the crystal 
plane. However, to experimentally observe the dependence 
of W with P [See Fig. 1(b)], we generated five more different 
states by changing the waist of the beam at the crystal plane. 
This has been done by using a configurable set of doublet 
achromatic lenses placed before the non-linear crystal [See 
Fig. 2(a)], The corresponding values of P for the six gener¬ 
ated states are shown in Tab. I. 


TABLE 1. Corresponding values of P for the generated states. 


State 

c [pm] 

P 

1 

200 

0.1595 

2 

100 

0.3189 

3 

70 

0.4556 

4 

45 

0.7087 

5 

40 

0.7973 

6 

35 

0.9112 


To guarantee the purity of the measured states, we used spa¬ 
tial and spectral filters in each measurement apparatus. For 
instance, interference filters were used to select degenerated 


In order to compute A 2 (<72 l^i), one shall perform condi¬ 
tional coincidence measurements at the far-held plane of the 
non-linear crystal. The reason is very simple: at this plane, 
the coincidence rate C q {x\ ,* 2 ) is given by 


C q {x\,x 2 )■ 


u 


■(*1 +X2) 


k 


■(*1 -X2) 


(9) 


and, since qj = kxj/f q , it maps the square of the probability 
amplitude of Eq. (7), i.e., the momentum correlation of the 
photons generated in the SPDC process. 

Figure 2(b) shows our setup configuration. One lens L q 
with focal distance f q = 15 cm was placed before the PBS to 
create the far-field plane for both signal and idler beams. At 
these planes, vertical slits were placed for post-selecting the 
desirable state [44, 45] and performing the conditional coin¬ 
cidence measurements [11, 31, 46]. The width of each slit is 
50 pm. The effect of using a non point-like detector is that 
the transmittance function of the slit may broad the far-field 
distribution to be measured. However, in the case of 50 pm 
slits and the experimental configuration adopted, this effect is 
negligible [42]. This can be easily checked through the calcu¬ 
lation of the convolution between the transmittance function 
of the slit and the predicted far-field distribution. 

After the transmission through each slit, the down- 
converted photons were collected with a lOx objective lens 
and multi-mode fibers. The fibers were connected to single¬ 
photon counting modules, and then a coincidence circuit (with 
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(a) 


(b) 



Position [m] |q 1 



Position [m] x 1Q -3 


(c) 


(d) 


1 



Position [m] x jq-5 



Position [m] x q 5 


FIG. 3. (Color Online) Experimental measurement of the far- and 
near-field conditional distributions. In (a) [(c)] we show the results 
for the second state at the far-field (near-field) plane. In (b) [(d)] we 
show the results for the sixth state at the far-field (near-field) plane. 
The (red) dots represent the experimental data and the dotted (black) 
lines are the theoretical predictions from Eq. (9) and Eq. (10). 


4 ns of coincidence window) recorded the data. To perform 
the coincidence conditional measurements of A 2 (q 2 \qi), we 
scanned in the horizontal direction one slit (of mode 2) while 
the other one was fixed at the center ( q\ = 0). 

To give an example of the results obtained while scanning 
the slit at mode-2, we show in Fig. 3(a) [Fig. 3(b)] the far-held 
conditional distribution measured for the second (sixth) state 
generated. The experimental results are represented by red 
points (error bars lie inside the points due to the observed high 
rate of coincidence counts) and the black dotted-line is the 
theoretical curve for these distributions arising from Eq. (9). 


B. Measurement of A 2 (x 2 |xi) 


To measure A 2 (x 2 |xi) it is now necessary to measure at the 
near-field plane of the non-linear crystal. Again, the reason 
is very simple: at this plane, the coincidence distribution is 
proportional to the square of the amplitude of Eq. (8), that is 


C x {xi,x 2 ) 



2 


( 10 ) 


For performing the conditional coincidence measurements 
at the near-field plane we used 4x objective lenses to form 
the image of the center of the BBO-II onto the transverse 
plane of the fiber-couplers [see Fig. (2)(c)]. For doing this, 
we removed the slits used to scan the coincidence rate in 
the far-field plane and the lens L c/ . The multi-mode fibers 
were replaced with single-mode fibers whose core diameter 
were 4.7 /. im. The small size of the fibers core allows for 
post-selecting and measuring with high accuracy A 2 (a' 2 |xi) 
[33, 40]. Again, the effect of the transmittance function of 
the fiber over the broadening of the near-field distribution is 



FIG. 4. (Color Online) (a) Experimental result and (b) theoretical 
prediction for the near-field conditional distribution while moving 
longitudinally the non-linear crystal (See the main text for details). 


negligible [42]. This can be easily checked through the cal¬ 
culation of the convolution between the transmittance func¬ 
tion of the fiber with the predicted near-field distribution for 
our experimental configurations. By scanning transversely 
the single-mode fiber at mode-2, we recorded the coincidence 
conditional distribution at the near-field plane. An example of 
the results obtained is shown in Fig. 3(c) [Fig. 3(d)] for the 
second (sixth) state generated. The experimental results are 
represented by red points. The black dotted-line is the theo¬ 
retical curve for these distributions arising from Eq. (8). 

We have measured the coincidence conditional distribution 
of the down-converted photons by imaging the center of the 
non-linear crystal to the transverse plane of the detection sys¬ 
tem. Flowever, as it has been shown in Ref. [33], the con¬ 
ditional distribution at the near-field plane depends strongly 
of which part of the crystal is imaged at the detection sys¬ 
tem. This is specially relevant when thicker crystals are con¬ 
sidered. In our case, we have checked that the near-field dis¬ 
tribution does not change significantly while considering dif¬ 
ferent planes of our thin crystal to be imaged at the detection 
system. For doing this, we moved the crystal around the lon¬ 
gitudinal central position zo = 0, imaging 5 different planes z c 
of the crystal at the detection plane. For each z c , we recorded 
the conditional distribution at the near-field plane. The exper¬ 
imental results are shown in Fig. (4). One can observe that 
our results are in agreement with the theoretical prediction, 
which takes into account our crystal length and our imaging 
system. A longitudinal crystal displacement around its center 
introduces a phase factor of exp[— i\cq\z/(2k)\ onto Eq. (7). 
From our results, one can see that there is only a slight nar¬ 
rowing of the near-field conditional distribution such that this 
effect does not affect significantly our test of W. 


C. Testing W 

In order to test the W witness, and then certify the non- 
gaussian feature of the CV spatial state of SPDC, we compute 
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TABLE II. Results of the conditional variances measured at the near- and far-field planes. 


P 


A 2 (-C 2 10 )e [m 2 ] 


A 2 fe|0 )e [ 


1 

wL 


W E 


Wi¬ 


0.1595 

(6.62 ±0.26) x 10- 10 

(3.55 ±0.14) x 10 7 

0.024 ±0.002 

0.033 

0.3189 

(9.17T0.37) x 10~ ll) 

(1.25 ±0.05) x 10 s 

0.115 ±0.009 

0.11 

0.4556 

(7.76±0.31) x 10~ 10 

(2.60±0.1) x 10 8 

0.20 ±0.02 

0.19 

0.7087 

(7.75±0.31) x 10~ 10 

(4.67 ±0.19) x 10 8 

0.36 ±0.03 

0.34 

0.7973 

(7.23 ±0.29) x 10~ 10 

(5.28±0.21) x 10 8 

0.38 ±0.03 

0.38 

0.9112 

(8.16±0.33) x IQ- 10 

(5.18 ±0.21) x 10 8 

0.42 ±0.03 

0.43 


the variances of the conditional coincidence measurements at 
the near- and far-held planes for the six generated states. The 
results of all the conditional variances measured are shown 
in Tab. II. The errors of the variances were obtained by min¬ 
imizing the squared two-norm of the residuals between the 
analytical and experimental results (See Fig. 3). 



P 


FIG. 5. (Color Online) Measurement of W for the six generated 
states. The continuous (black) line corresponds to the theoretical 
value of the inequality for the spatial SPDC state. The dashed (blue) 
line shows the Gaussian bound. 

Figure 5 shows the values of the W values for each state 
generated. The continuous (black) line indicate the predicted 
theoretical value of the inequality when we consider the spa¬ 
tial state of SPDC. The dotted (blue) line is the upper limit for 
this inequality for Gaussian states. For the first three entan¬ 


gled states, the values obtained for this inequality are below 
the Gaussian limit as predicted by theory. However, for the 
last three entangled states, there is a clear experimental viola¬ 
tion of this bound. Due to momentum conservation, the spa¬ 
tial state of the SPDC process is always entangled [41] [See 
Fig. 1(a)] and, thus, we have a clear experimental demonstra¬ 
tion of the non-Gaussianity of the CV two-photon spatial state 
of SPDC. 


V. CONCLUSION 

We have introduced a novel application for the EPR- 
steering inequality by showing that it can be used for witness¬ 
ing the non-Gaussianity of CV quantum states. To demon¬ 
strated this we performed an experiment using the CV spatial 
state of entangled down-converted photons. Due to the phase¬ 
matching conditions of the SPDC process, the generated is 
state is naturally a pure entangled non-Gaussian state. A clear 
violation of the Gaussian bound of the EPR-steering inequal¬ 
ity has been observed. Since non-Gaussian states are required 
for many new protocols of CV-QI, our work highlights the 
simplicity and relevance of using SPDC sources for new ap¬ 
plications in CV quantum information processing. 
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